Effect of next-nearest neighbor coupling on the optical spectra in bilayer graphene 
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We investigate the dependence of the optical conductivity of bilayer graphene (BLG) on the 
intra- and inter-layer interactions using the most complete model to date. We show that the next 
nearest-neighbor intralayer coupling introduces new features in the low-energy spectrum that are 
highly sensitive to sample doping, changing significantly the "universal" conductance. Further, its 
interplay with interlayer couplings leads to an anisotropy in conductance in the ultraviolet range. 
We propose that experimental measurement of the optical conductivity of intrinsic and doped BLG 
will provide a good benchmark for the relative importance of intra- and inter-layer couplings at 
different doping levels. 
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PACS numbers: 73.50.Mx, 78.66.-w, 81.05.Uw 

Since the isolation of single layers of graphite in 2003 
1], a lot of exciting work on single layer graphene (SLG) 
has been done |2j . For example, the prediction and obser- 
vation of electron-hole symmetry and a half- integer quan- 
tum Hall effect [1,13, 0L_ finite conductivity at zero charge- 
carrier concentration [3j, the strong suppression of weak 
localization [1, 0, [1], universal conductance [H, [l(| U| 
and magnetic enhancement of optical conductance in 
graphene nanoribbons [12j |. 

More recently, attention has also been paid to SLG's 
cousin, bilayer graphene (BLG). The electronic and 
transport properties of BLG differ significantly from SLG 
in many respects, particularly at low energies in the 
'Dirac' regime. Various models for low energy BLG ex- 
ist in the literature depending on the coupling terms in- 
cluded, and whether electronic bands beyond the lowest 
energy subbands are retained [3, [3|. Many interest- 
ing results were obtained based on a model that includes 
only the most dominant of the interlayer coupling terms 
in BLG, as well as the usual nearest neighbor intralayer 
term 15]. By including the second most dominant inter- 
layer coupling, some unusual properties such as a peculiar 
Landau-level spectrum have been derived [lij , as well as 
a new low energy peak [l(| . 

The 'universal conductance' of graphene is both a DC 
and an AC phenomenon. It is a direct result of the lin- 
ear energy dispersion of graphene. Linear subbands im- 
ply both a constant density of states as well as consis- 
tent transition matrix elements, which means that for 
as long as the linear (Dirac) approximation is valid, the 
conductance is a constant. In the AC case, the value 
of the universal conductance of single layer graphene is 
a i = e 2 /Ah. In the layered standard benchmark is 

simply o~ n = no\. However, this is not generally accurate, 
as the subband curvature caused by interlayer coupling 
in the case of layered graphene leads to a non-constant 
conductivity. This raises an important question: in what 
energy range is a n = no\ applicable? 

The infrared conductance of BLG has been measured 



by several groups [17|, [18( . These results rely upon the 



effects of an induced gate voltage on the bandstructure, 
and all assume a discrepancy in onsite energy between the 
two layers. In ref. 17, Mak et al present the 'expected' 
IR conductance without the latter assumption, and find 
that it differs markedly from their experimental results. 
This demonstrates the need to assume an energetic dis- 
crepancy between the two layers in BLG. Our theoretical 
results, however, show a strong correlation to the results 
in ref. 16 and 17, demonstrating that while an energetic 
discrepancy may exist, it is not necessary in describing 
the IR response observed experimentally. 

In this letter, we study the dependence of the optical 
conductance of BLG on various intra- and inter-layer cou- 
plings. It is shown that the interplay of these couplings 
leads to a significant deviation in the behaviour of the 
conductance at low frequencies, which can, in turn, be 
tuned by electronic doping. In the important ultraviolet 
frequency band, this interplay leads to significant con- 
ductance anisotropy, i.e., the absorption along the arm- 
chair direction is around 50% stronger than that along 
the zigzag direction. 

A typical BLG sheet consists of two SLG layers stacked 
in the orientation shown in Fig.l. Several forms of the 
Hamiltonian for BLG are used in the literature depending 
on the approximations used and the relative orientations 
of the two layers. The original consideration was given 
by Slonczewski-Weiss-McClure which included all three 
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interlayer coupling terms 



The most prominent 



interlayer term is the A-B and B-A coupling between sites 
which are directly above (or below) each other. Here we 
define this term as 71 = 0.36ey. The other two inter- 
layer coupling terms are the A-B and B-A coupling be- 
tween inequivalent sites which are not directly above or 
below each other, but offset by an amount b — 1A2A, 
and the A-A and B-B terms which are similarly offset 
from one another, but represent equivalent sites in the 
SLG Brillouin Zone. These coupling terms are defined 
here as 73 = OAOeV and 74 = 0A2eV, respectively. We 



FIG. 1: The three interlayer and two intralayer coupling terms 
included in the BLG Hamiltonian. 73 and 74 differ in that 
they connect, respectively, inequivalent (eg. A-B) and equiv- 
alent (eg. A-A) points in the SLG Brillouin Zone. 71 is a 
directly vertical transition, and so the overlap of the wave- 
functions is about 3x larger than 73 and 74. The armchair 
(AC) direction is given by the x-axis, and the zig-zag (ZZ) 
direction is given by the y-axis. The lattice vectors ai and a.2 
are also shown. 



have also included the next nearest neighbor A-A and 
B-B coupling which we define as t' = 0.30eV. Finally, 
as usual, the nearest neighbor A-B and B-A coupling is 
included, which is given here by t — 3.0eV . All energies 
will be normalised relative to the first nearest neighbor 
coupling. 

The full Hamiltonian matrix for the BLG system is 



H 



BLG 



( t'H' tH* 7 4 i? 71 \ 

tH t'H' -y 3 H* 74-ff 

7 4 ff* j 3 H t'H' tH* 

\ 71 74-ff* tH t'H' J 



(1) 



where H = e^Wv^i + e * a+ + e ik-a_) and jji = 
2(cos(k • a + ) + cos(k • a_) + cos(k • (a+ — a_))). Here 
a± = a(±|,— ^) are the two lattice vectors shown in 
Fig.l. The eigenvalues and eigenvectors in the absence 
of 74 are readily solved. With 74 included however, the 
form of the solution is unwieldy. The eigenvalues in the 
simpler case are given by the (relatively) concise form 



i'(esLG - 3) + a 




(2) 



Where 
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FIG. 2: The k x dependence of the bandstructure near the 
K/K' points. The red dashed arrow represents a transition 
which is permitted in an undoped sample if NNN coupling is 
neglected, but becomes forbidden when it is. The black solid 
arrow is the opposite: a previously forbidden transition be- 
comes allowed when NNNs are included. The effect of doping 
is to raise or lower the Fermi level, making the inclusion of 
NNNs partly equivalent to doping. 



And 7j^ = 7fe|j, ± 7i, with s,s' — ±1, and esL are 
the regular eigenvalues for the SLG system given as 



£ slg = t(l + 4cos(afca;/2) cos(ak y /2V3) 
+ 4 cos 2 (ak x / '3)) \ 



(4) 



(3) 



From this result we see that there are two conduction 
bands and two valence bands which are confined above 
and below the line e s . s ' — £'(egLG — 3), respectively. This 
simple result will form the basis for much of the discus- 
sion to follow. 

The electron field operators can be con- 
structed from the eigenvectors ip StS i(]c) such that 
tf(r) = (l/47r 2 )E k , s , s ,a s ,,'(k)^j'(k)e lk - r , where 

a SjS '(k) (a s s i (k)) denotes the annihilation (creation) 
operator for an electron in the s or s' subband with 
momentum k. 

The band structure of BLG near the K points varies 
dramatically depending on the coupling terms included 
in the Hamiltonian. The effect of the various coupling 
terms are as follows: 

The next nearest neighbor term t' defines the line 
about which the conduction and valence bands are con- 
fined, as can be seen in Fig. 2. In the absence of lateral 
warping, the valence and conduction bands are symmet- 
ric about this line. Near the K points, the prefactor of 
the energy term associated with t' approaches 3, and so 
the effect of this term is to shift all four bands downward 
by an amount 3i' w O.leV. This can easily be seen from 
equation 2, where H' w 3. This breaks the electron-hole 
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FIG. 3: The k x dependence of the two inner bands near the 
K/K' points zoomed right in to see the effects of the intralayer 
next nearest neighbors (NNN), and interlayer coupling terms 
73 and 74. The NNN interaction has shifted these features 
well below the Fermi level. 73 causes a second dirac point to 
emerge, and 74 skews the bandstructure, causing one of the 
two Dirac points to be pushed down to a lower energy. 



symmetry, but not the location of the minima in momen- 
tum space. 

The term 71 represents the dominant interlayer A — B 
and B — A coupling. This term causes an energy gap to 
form between the two conduction bands, and an identical 
gap between the two valence bands of 71 ss 0.13eV '. 71 
also removes the linear dispersion at low energies. The 
electron hole symmetry is retained, and no lateral warp- 
ing occurs. The effect of 71 is apparent in Fig. 2. 

The second interlayer coupling term 73 restores the 
linear lowest energy subband. This term causes what 
is usually referred to as 'trigonal warping' [14|, |22|. A 
second set of Dirac points near the K/K' points emerges 
with 73 included, as can be seen in Fig. 3. 

Finally, 74 causes one of the Dirac points to be plunged 
below the NNN line, also seen in Fig. 3. When the next 
nearest neighbor term couples with 74, however, the low 
energy x-y isotropy is substantially weakened. While 
73 causes the well known anisotropic 'trigonal warp- 
ing', the energy range of this effect is in the order of 
i/10000 0.0003eF. On its own, the effect of 74 is simi- 
larly small. Here, however, 74 and t both couple equiva- 
lent sites, which causes a compounding of their individual 
effects on the electronic dispersion relation. The effect is 
quite a large deviation from isotropy. The effect of this 
deviation is most noticeable in the low energy conduc- 
tance anisotropy shown in Fig. 5. 

We now evaluate the optical conductivity of BLG in 
the absence of disorder or impurities, over all relevant 
photon energies. By using the Kubo formula, the optical 
conductivity is given as [23j, 
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FIG. 4: The full energy optical conductance (in units of 
<7i = e 2 /h) vs the normalized frequency Q = foui/t for bi- 
layer graphene. Generally, a xx (the armchair direction) has a 
larger optical response than a yy (the zig-zag direction). When 
NNN and 74 are neglected, and at low energies, a xx = a yy . 
This is no longer the case here, with NNNs and 74 included. 
For e < 3t' , there are no allowed transitions and the OC is 
approximately zero. The grey shaded area indicates the low 
energy region plotted in Fig. 5. 



^H = -/ dte^dJ^MO)]). (5) 
w Jo 

The components of the current operator can be calcu- 
lated from J v ^{t) = e tHt J^(0)e- tHt , where J^(0) = 
* t (r)w„, A1 v E'(r'), in which % jfl = dH/dk u ,n, and v, fi = 
x,y. These values are calculated numerically, but we 
note that for each band there are three types of inter- 
band transitions, and also intraband transitions. In the 
case of no disorders and no intermediate interactions, it 
is found that intraband transitions cannot occur. 

In Fig. 4. we examine the full energy optical conduc- 
tance of undoped bilayer graphene with NNN interac- 
tions included. Near the higher energy valley points, the 
optical conductivities exhibit two extrema, similar to the 
single peak found in single layer graphene [2l|, [24[ • These 
peaks correspond to the two dominant vertical transitions 
between the two symmetric pairs of saddle points. The 
Joint Density of States in these valleys reaches a cusp-like 
maximum which leads to the extrema in the conductivi- 
ties. These two energy peaks are separated by an amount 
hw = 271, as expected from the bandstructure calcula- 
tions. 

Fig. 5. shows the low energy optical conductivity of the 
undoped sample from the grey shaded region of Fig. 4., 
as well as a sample doped to the level of the NNN in- 
teraction. The longitudinal conductivities vary greatly 
when including a non-zero NNN interaction. The ef- 
fect of the dominant interlayer term 71 (ie. setting 
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FIG. 5: The low energy optical conductance region from 
Fig. 4. at two different doping levels. The black solid line is 
the OC of undoped BLG. The blue dotted line represents the 
OC at finite doping, and is equivalent to an undoped sample 
with no NNN interaction included which has been previously 
reported. The NNN-74 coupling causes a new peak to emerge, 
and suppresses the previously reported one. This new peak 
is much larger and shifted to a lower photon energy. In a 
suitably doped sample, however, the t' = (no NNN) peak 
has been retrieved by an effective shifting of the Fermi level. 



73 = 74 = t' = 0) at low energies has been reported 
recently This result has been retrieved in our re- 

sult by doping the sample to /1 = —O.leV/t. In the 
undoped curve however, the previous result is entirely 
suppressed, and replaced by an approximately 2x larger, 
rounded peak, followed by a significant trough. This cor- 
relates well with the behaviour observed in IR experi- 
ments 



171 ll8[ , although without the added effects of an 
induced gate voltage. The doped case represents a transi- 
tion of electrons from the upper conduction band to both 
the upper and lower valence bands, which is equivalent to 
neglecting the NNN coupling, as shown in Fig. 2. The un- 
doped case however reflects the suppression of transitions 
into the lower valence band, since it is already filled, and 
yet a new set of transitions occur between the two valence 
bands. These bands are separated by an approximately 
constant factor of 71 , which leads to the large peak cen- 
tred at Jloj ~ 71 in the undoped bilayer case. The feature 
is in striking contrast with that of SLG. For SLG, the 
effect of the NNN coupling is to suppress the universal 
conductance at low frequencies 24 1 . As is clearly seen in 
Fig. 5., for BLG the interplay of the interlayer coupling 
and the NNN coupling can suppress the conductance at 
low frequencies. However, it also induces a strong ab- 
sorption peak in the far infrared before the onset of the 
universal conductance. 

For this reason, the low energy approximations of the 
behaviour of bilayer graphene are generally more relevant 



to doped samples, with the undoped bilayer properties 
being drastically affected by the next nearest neighbor 
hopping and additional interlayer terms. Furthermore, 
when using existing theories to explain experimental re- 
sults, it needs to be noted that an energetic discrepancy 
between layers, as well as the inclusion of a gate voltage, 
both cause some similar effects to the inclusion of the 
NNN interaction. All of these will therefore need to be 
accounted for when explaining any experimental result. 

Furthermore, the conductance anisotropy observed in 
Figures 4 and 5, which is prominent even in the IR region 
when 74 and t' are both included, makes the polarization 
of the photon beam in experiments a relevant parameter. 
This orientation dependence of the optical conductance 
makes determination of the orientation of a BLG flake 
possible, and also makes BLG a potential partial polar- 
izer. The doping dependence of the low energy conduc- 
tance anisotropy makes this feature quite versatile. 

Finally, as we have already mentioned, the value of 
the 'universal' conductivity is a topic of great interest 
at the moment. According to these results, which have 
been calculated from the most robust interlayer and in- 
tralayer model adopted to date, the value of the universal 
conductivity is a = 2u\ where o\ = e 2 /Ah is the univer- 
sal optical conductivity of single layer graphene defined 
earlier. The range over which this value is applicable 
is greatly affected by the presence of the NNN interac- 
tion, and the electronic doping. In particular the NNN 
interaction causes the very low energy optical response 
to become negligible, and around the observed peak, the 
optical conductance is <T poa k ~ 4ai. 

In conclusion, we have studied the longitudinal opti- 
cal conductivity of BLG with the inclusion of all relevant 
interlayer coupling terms and next nearest neighbor in- 
tralayer interactions. The optical conductivity exhibits 
double peak resonance separated by an amount 271 and 
centered around tiw = 2t. At low energies, the NNN 
interaction leads to entirely new behaviour of the op- 
tical conductivity. The results obtained without NNN 
coupling, however, can be retrieved by appropriate elec- 
tronic doping. The interplay of the NNN-74 couplings 
were found to lead to significant low energy conductance 
anisotropy which is strongly doping dependent. Finally, 
the value of the universal conductivity with the most ro- 
bust formalism used to date has also been determined, 
and is given by 02 = 2<7i . These results will be crucial to 
the experimental testing of accepted theories on bilayer 
graphene, and will be useful for potential low energy elec- 
tronic and photonic applications of bilayer graphene. 
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